Effective theory of fluctuating orbital currents in high-T c cuprates 
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We derive an effective dissipative quantum fieid theory for fluctuating orbitai currents in clean 
C11O2 sheets of high-T c cuprates, based on a three-band model. The Coulomb repulsion term 
between Cu- and O-sites is decoupled in terms of current operators representing horizontal and 
vertical parts of circulating currents within each C1LO2 unit cell of the lattice. The model has 
ordering of currents at finite temperatures. The dissipative kernel in the model is of the form 
|w|/|q|, indicating Landau damping. Applications of the effective theory to other models are also 
discussed. 
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Constructing an effective description of the long- 
wavelength and low-energy physics of high-T c supercon- 
ducting cuprates represents a profound and formidable 
problem in physics. Such a description must be consis- 
tent with experimental observations of several anoma- 
lous normal state properties of these systems. Varma 
has recently proposed that quantum critical fluctuations 
associated with the breakup of a subtle order, involving 
circulating currents, could induce the observed anoma- 
lous normal state properties of high-T c superconductors 
Essentially, the associated quantum critical fluctu- 
ations are suggested to produce a fluctuation spectrum 
resulting in a Marginal Fermi Liquid [2]. Recently, such 
a spectrum has been derived from a conjectured effec- 
tive field theory of circulating currents [3j. It should, 
however, be mentioned that a recent numerical evalua- 
tion of the current-current correlations in a three-band 
t — J-model with 24 sites, where doubly occupied sites 
have been projected out, shows no evidence of the orbital 
current pattern proposed by Varma [4j . 

The particular form of proposed order involves circu- 
lating currents within a CuOi unit cell where the cur- 
rents run horizontally and vertically through a Cu site 
and close by direct hopping between O orbitals, as in Fig. 
[TJ Three other equivalent patterns may be found by re- 
versing the direction of the current through each Cu-site 
in the horizontal and vertical directions. This results in 
a pattern of staggered orbital magnetic moments within 
each unit cell, such that the pattern repeats from unit 
cell to unit cell. A magnetic intensity of the type associ- 
ated with the above orbital magnetism has recently been 
reported [f|. Since no obvious thermodynamic singular- 
ities have so far been reported at the pseudogap-line in 
the cuprates, it is important to investigate whether or 
not the proposed models for this novel type of order im- 
ply the presence or absence of prominent signals in such 
quantities as specific heat or (indirectly) magnetization. 
Other staggered orbital magnetic patterns have also been 
proposed, most notably the extension of the staggered 
flux phase [|| to finite doping 0] ■ We emphasize that our 
justification for focusing on the orbital current pattern 




FIG. 1: (Color online) The circulating current phase Qu pj|. 
Cu sites are grey circles, O sites are red. The unit cell is 
shown by the dashed square. A staggered magnetic moment 
pattern within each unit cell that repeats from unit cell to 
unit cell (the curl of the blue directed circles) is indicated. 



proposed by Varma [![, are the experiments reported in 
Ref. i. 

We derive an effective quantum field theory for fluc- 
tuating orbital currents from a microscopic descrip- 
tion of clean CuOi planes. We are primarily inter- 
ested in investigating the intrinsic effects such fluc- 
tuations have on the physics of the cuprates. We 
therefore neglect disorder, as was also done in Ref. 
3. With ever improving sample quality, we expect 
that the effective theory we derive should be use- 
ful. The starting point is the three-band model H = 



H 



(i) 



H^l where 



Y,r,a £ d d l,a d r,a + K pd + K pp , 
t^pd V Sr,<j Nr,ff (?i,r+f x,<r Px,r— §x,ct Py,r+^y,a + 
Py,r-§y,<?) + h.C.], K pp — — t pp X}r,a-[(P:c,r+§ &,a ~~ 

i , !,r-§*,«r)(Pv.r+t*.* ~ Pv,r-&,<r) + h ' C '] and H int = 
^ X/r,CT,cr' n d,r,cr{ n p x ,r+§&,a-' 4" 1p„r- -fx,<7' + n p y ,r+ §y ,a' + 

n p y ,r--§ y.o-')- We work with electron operators and the 
vacuum is defined as empty d x i_ y i, p x and p y orbitals. 
The r-sum runs over the Cu-lattice. The Cu-O and O- 
O hopping is governed by the parameters t pc i and t pp , 
respectively, whereas e e i is the difference in on-site en- 
ergy between the copper and oxygen orbitals. The term 



2 



i? int represents on-site repulsion terms, for which we 
make the crude assumption that their effect is to merely 
renormalize the hopping parameters t pc i — > i pc t = \x\t p d, 
t pp — > t pp = \x\tpp, where is the deviation from 
half-filling [l|. We also assume the O-O repulsion to 

(2) 

be small. Hence, we only consider explicitly H int ', the 
Cu-O-repulsion. 

The interaction-term can be decoupled in 

terms of bosonic fields coupling to the bilinear fermion 
operators A« >(T , = £ k «L- q pl ,k- q ,a> + 

a y,k-q^,k-q, CT ') d k, CT with i = 1,..,4 . Here, N is the 
number of Cu lattice sites. We define a^\. — = 
siri(k x a/2) = 8 X) k, <x^ k = = cos(fc K a/2) = c x<k and 

(1) (2) ' (3) (4) , 

a U = - a v ,k = s v,k, %,k = ~ a y,L = c v,k, where a is 
the Cu-Cu lattice constant. A discussion of {A^\ a ,) as 
translational invariant order parameters in the cuprates 
is found in While i = 2 transforms as the kinetic 
energy, i — 1 and i = 3, 4 give rise to different current 
patterns. Since the observed magnetic signal [5| is con- 
sistent with the current patterns of i = 3, 4, we keep 
only this in what follows. An effective model for the 
i = 1-part was considered in Q. Observe the relation 

where, in real space, 

3r,<T,*> = 2 [ d r,° (^r+f*,^ + fc,r-f x,a') ~ h.C.] . (1) 

In a unit cell centered on Cu, this is proportional to the 



current from the left oxygen to the copper plus the cur- 
rent from copper to the right oxygen. We define j v v a a , in 
the same way, but with a minus sign due to the d-wave 
symmetry of the Cu-orbital. Finite expectation values of 
«r ct <t' would correspond to Landau-Pomeranchuk insta- 
bilites, believed not to be relevant in the cuprates. Thus, 
we retain only the decoupling fields that correspond 
to spin diagonal expectation values of the operators 
j* ^/(t), since (jr^aij)) 7^ in the current pattern 
depicted in Fig. [TJ The fields retained, Jt^ v \t), are 

real and (Jr (y) (r)> = ^(ir,? CT '( r )><W'> «- e - thc ficlds 
represent charge currents on horizontal and vertical 

O-Cu-O-links. The fields J* (t) and the fermions 
are coupled by particle-hole excitations of thc form 

^Ek.q^ (yv^-^i.k-q.^A^ - o -> y) - h.c.), 

where the time dependence was omitted. It is important 
to keep in mind that the bosonic fields Jr (t) transform 
as vectors under a change of coordinate system. Note 
that we could also have chosen the arguments of the 
a^'s to be k and not k — q in a , , corresponding to 
a decoupling in terms of currents defined on horizontal 
and vertical Cu-O-Cu-Ymks. 

Integrating out the fermion fields, we obtain 
the partition function as Z = JDJ X DJ V e~ s , 
where the effective action is given by S = 

W £q,0,„ [Jq(^)^q(-W !/ ) + J*(iw v )JH.q(-iU v )] ~ 

Tr In [Gq 1 + S] . Using the gauge transformation 
PxM,<y -> iPx,h.,<r, Py.Kcr -> -iPy,\t,<T> we have 
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where we have defined ki2 = ki — k2 andwi2 = uj ni —uj n2 . 
For tpp = 0, the non- interacting part of the problem Qq 1 
may easily be diagonalized into three quasiparticlc bands 
4 0) = 0, 4 ±} = e„/2 ± V / ( £d /2)2 + 4^( s 2 ifc + s 2 fc)i of 

which # k 0) , E^ 5 are full and E { k +) is partially filled. This 
picture is not qualitatively altered by t pp ^ 0. A nonzero 
value of t pp is however vital for the realization of the 

current pattern. It is implicit that (J? (t)) — > when 
Expanding the last term [lfj| . odd powers of J van- 



ish, such that Trln [G' 1 + E] = Trlng" 1 - ±Tr [£ £] 2 + 
C(J 4 ), where — TrlnC?^ 1 gives the free energy of the non- 
interacting system, and £ involves the fluctuating fields 
J x and J v . To second order in the fields J X W and in 
space and imaginary time gradients, we have derived a 
quantum dissipative effective action S = Sc + Sq, where 

Sc = GclijJqiiu^Jl^-iuu), 

Sq = E E G Q^ J q(^) J -q(-^)- ( 4 ) 

q,Wt, i,j=x,y 
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with G C L XX = a c + ai q 2 x +a t g 2 , G c ^ yy = a c +a h q y 



G 



C.xy 



G 



C,yx 



> 



xyQxqy, G Q>XX = a wl+ a d ^ql, 

^-WVl aild G Q l xy = ^yx 

Here, q x = q x /\q\- The dissipation ker- 
nel is valid for |w v |/|q| <C 1. The limit |w„|/|q| 3> 1 does 
not contribute to dissipation. The explicit expressions 
for the coefficients cti are unwieldy and of limited use. 
The equality of the diagonal and off-diagonal dissipation 
coefficients is only correct when t pp = 0. Changes when 
t pp =/= are small and unimportant, and are neglected 
in the following. Note also that this theory might not 
be applicable to the ordered phase, since the Fermi sur- 
face is proposed to be gapped there However, it is 
the fluctuation spectrum in the disordered phase which is 
important in connection with the Marginal Fermi Liquid 
hypothesis [2J. 

We have divided the action into a classical (C) and 
a quantum (Q) part. At finite temperatures, only the 
classical piece of the action Sc needs to be considered as 
far as critical properties are concerned. The excitation 



energies of the eigenmodes of Sc are given by X± = a c 
(ai +a t )q 2 /2± yJ(Aa) 2 q 4 + "iq 2 x q 2 , where Aa = (a; 
= a t )/2, and 7 = a 2 - (2Aa) 2 . Hence, for (a/,a t ) > 0, 



v 

uniformly ordered state is stable in the classical domain 
below some critical temperature, provided a xy < af+a 2 . 

The dissipation kernel essentially gives Landau damp- 
ing, albeit anisotropic due to the directional nature of the 
fields. The dissipation is a result of coupling to the gap- 
less particle-hole excitations in the band -E^ • The sm ~ 
gular form |av|/|q| is correct only if the order in the hor- 
izontal and vertical currents are uniform and not modu- 
lated at some nonzero reciprocal vector. It implies that 
the dynamical critical exponent z = 3 12|. See however 
Ref. M 

Current amplitude fluctuations are expected to be 
high-energy excitations [l[ and will therefore not deter- 
mine the critical properties of the model. Thus, we treat 
the fields J, (r) and J%{t) as Ising variables. Reverting 
to a real space Cu-lattice formulation and setting a = 1, 
we obtain (up to constant terms) 



Sc 



Sq = a a 



ch 



(& X r,r'J:(r)J^(T)+al ir jy(r)J y r ,(r))+ £ fi£, ( J* (r) J* (r) + J r » J$ (7 



<r,r'> 



«r,r'» 



(5) 
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gjg \ 2 fdjy\ 2 



dr J \ dr J 



r r' i 4 



Here, (r, r') and ((r, r')) denote nearest-neighbor and 
next-nearest-neighbor summations, respectively. For 
r — r' = ±x, a^ r , = &i and a y r , — at, whereas 
when r — r' = ±y, a* r , = at and a y rr , — &i. The 
parameter a^ v r , — a xy when r — r' = ±(x + y) and 
a* y r , = — ct xy when r — r' = ±(x — y). The coef- 
ficient ad > and the positive semidefinite matrix 
K r _ r '(r - t') = K r - r *(T - r')g r _ r / ® g r -r', where 
g r _ r , = (r - r')/|r - r'| and K r {r) = l/(|r| sin 2 (7rr//3)). 
Fluctuations ( —J^,Jf — > J%) corresponds to going 
from the depicted current pattern (Fig. [1]) to a new one 
which is obtained by a counterclockwise rotation by ir/2, 
(J, — > Jr,J$ — * —J y ) corresponds to clockwise rotation 
of 7r/2, and (J* — > — , — > — J^) to a rotation of 7r. 
It is implied that in the dissipation kernel, we must use 
a short-distance cutoff in (t, r)-space, since the expres- 
sions are derived in the limit of low (w,q). The tildes 
on the coefficients indicate that the model in Eq. [5] is 
regularised on a lattice, and that the fields have been 
normalised to Ising- variables. Moreover, there will be 
higher order (quartic) terms generated that simply in- 
volve local squares of Ising variables multiplied by som 



bilinear combinatoin of Ising variables, and these will also 
contribute to the coefficients of the quadratic terms even 
before a renormalization group analysis is carried out. 
These terms are also taken into account by the tilde. 

In general, we have &i ^ ctt- A current living on a 
horizontal O-Cu-O-link, J*, couples to Jr±x through at, 
and to J£±y through at ■ As seen from Figure [fl there is 
no reason for these couplings to be similar, and in fact a 
detailed derivation shows that they are not 

At finite temperature, we may ignore the inertial and 
dissipative terms, which reduces the model to a classical 
model of two coupled Ising fields. Such a classical model 
will suffice to study the breakup of the current pattern 
at finite temperatures, while its quantum critical version 
can only be accessed via the full dissipative field theory. 
Note also that the dissipation kernel is non-local both in 
imaginary time and in space. The latter distinguishes 
this dissipation term from the Caldeira-Leggett type of 
dissipation appropriate for an array of Josephson junc- 
tions [13L 1 1 51 ] . The non-locality in r-space is anisotropic 
for the same reason as for the nearest-neighbor coupling. 

Eq. ([5]) may be rewritten on the form 
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Sc 



2ao / dr 
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r,r' l,j 
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/ dTdT 'EE ( Jl ( r,r) - ^(^,rO)K^ r ,(T- r')(j J '(0 r ,r) - J J '(^,r')) . 

^ V I*' 7 l" 



where we have used the parametrization cos(# rr ) = 
(J?(r) + J?(r))/2, sin(^ r ) = (J*(r) - J y r {r))/^ and 
#r,r € (0, 7r/2, 7r, 37r/2). We have denned a = (a; + at), 
(Aa) r ^' = (a; — &t) for r — r' = ±x and (Aa) rir / = 
-(dtj - a t ) for r r' = ±y. 

Eqs. ([5]) and © are the main results of this paper. 
These models describe a phase transition from a disor- 
dered bosonic state (a Fermi liquid), into a state with 
bosonic order in the form of ordered orbital currents. 

We next proceed to discuss some qualitative aspects. 
Consider first this model at finite temperature, where we 
may use the approximation S ~ Sc- When (ai,a t ) > 
and a xy = 0, the current pattern in Fig. [T] repeats uni- 
formly from unit cell to unit cell throughout the system 
in the ordered state. The specific heat has a logarithmic 
singularity at a critical temperature determined by the 
condition sinh(2/? c aj) sinh(2/3 c at) = 1, where f3 c — l/T c . 
Anisotropy in the nearest-neighbor couplings suppresses 
the critical temperature and critical amplitudes, and nar- 
rows the critical region, but does not alter the universal- 
ity class of the phase transition When (a;, at) = 
and a xy ^ 0, the ground state of the system features a 
striped phase in the diagonal directions, irrespective of 
the sign of a xy . Note also that when sign(a ; ) ^ sign(d t ) 
and a 



xy 



0, one obtains order with a period of twice 
the lattice constant. 

The dissipative term in this model comes from the cou- 
pling of the bosonic current fields to particle-hole exci- 



(+) 



i.e. an intra- 



tations in the partially filled band Ey 
band transition. In the above, we defined the currents 
on horizontal and vertical O-Cu-O-links, living on Cu- 
sites. We could alternatively have defined the currents on 
Cu-Cit-links, both in a three-band model and in a one- 
band model. This definition would be relevant to the 
study of d-density waves 0, 0]- However, one would ex- 
pect a different dissipation term in that case, due to the 
finite modulation vector of the ordered currents. Note 
also that the w^-terms in Sq in Eq. (U|), equivalently 
the inertial terms in Eqs. (O and ©, are of multiband 
origin. 

The quartic terms in 5c that would emerge from the 
above treatment are of the type cx l J^ T3T4 J^ 1 J} 2 Jl 3 J^. 



ijlm 
rir 2 r 

Note that for i = j = x, I = m = y, these terms in- 
clude an Ashkin- Teller type of four-spin interaction, used 
in Ref. y to argue that the Ising type of singularity in 



specific heat would be quenched. 5c in Eqs. (|5|6p dif- 
fers from the model of Ref. in several respects. How- 
ever, a direct comparison is difficult, as it is not clear 
what physical quantitites the fields in Ref. represent. 
Firstly, the Ising-exchange coupling terms in Eqs. (|5|6j) 
are anisotropic, possibly highly anisotropic, due to the 
bond-character of the Ising variables. Moreover, the term 

& r V r> J ? J r> in E( 1 S - <GH> is abseilt in Ref ' S While 

this term may be perturbatively irrelevant, it is far from 
clear that a xy is actually small. In addition, there also 
seems to be a discrepancy between the dissipation kernel 
|w„|/|q| derived here and the one employed in Ref. S 

We expect our model to be generically useful in de- 
scribing thermal and quantum critical fluctuations of di- 
rected particle-hole bond variables in fermionic lattice 
models. 
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